Screened QED corrections in lithiumlike heavy ions in the presence of magnetic fields 
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A rigorous evaluation of the complete gauge-invariant set of the screened one-loop QED corrections to the 
hyperfine structure and g factor in lithiumlike heavy ions is presented. The calculations are performed in both 
Feynman and Coulomb gauges for the virtual photon mediating the interelectronic interaction. As a result, the 
most accurate theoretical predictions for the specific difference between the hyperfine splitting values of H- and 
Li-like Bi ions as well as for the g factor of Li-like Pb ion are obtained. 



PACS numbers: 31.30.jf, 31.30.Gs, 31.30.js 

Investigations of the hyperfine splitting and the g factor in 
highly charged ions give an access to a test of bound-state 
QED in strongest electromagnetic fields available at present 
for experimental study. To date, accurate measurements of 
the ground-state hyperfine structure and of the g factor were 
performed in H-like ^o^Bi, i^Sho, ^^^Re, i^^Re, 207p^^ 203-^^ 
and 205T1 and in H-like [g] and ^^O H, re- 

spectively. In particular, the 2002 COD ATA value for the elec- 
tron mass is derived mainly from the experimental and the- 
oretical g factor values for hydrogenlike carbon and oxygen 
with an accuracy 4 times better than that of the 1998 CODATA 
value. An extension of such kind of experiments to highly 
charged Li-like ions presently being prepared will pro- 
vide the possibility to investigate a specific difference between 
the corresponding values of H- and Li-like ions, where the 
uncertainty due to the nuclear effects can be substantially re- 
duced iflOi [Tl , 12 1 . Achievement of the required theoretical 



accuracy for the hyperfine structure and for the g factor in the 
case of Li-Uke ions is a very interesting and demanding chal- 
lenge for theory. 

At present, the theoretical accuracy of the specific dif- 
ference of the hyperfine splitting values of H- and Li-like 
ions and of the g factor of Li-like heavy ions is mainly lim- 
ited by uncertainties of the screened QED and higher-order 
interelectronic-interaction corrections. In the present Letter 
we focus on one of the most difficult correction, namely, 
the screened QED correction in the presence of a magnetic 
field perturbation. State-of-the-art evaluations of the screened 
QED correction were performed with local screening poten- 
tials [IT3, 14, isj]. These calculations are based on the well- 
established technique developed for the evaluation of the one- 
loop QED coiTections in the presence of an external potential 
III6II . However, the employment of a local screening poten- 
tial does not allow one to take into account consistently all the 
contributing diagrams and to provide a reliable estimation of 
the uncertainty of the result. Therefore, a systematic descrip- 
tion in the framework of QED requires the use of perturbation 
theory. This crucial step has been made now and in this Let- 
ter we report on our results of the rigorous evaluation of the 



complete gauge-invariant set of the screened one-loop QED 
corrections. As the most interesting application of these re- 
sults towards tests of the magnetic sector of bound-state QED 
we present improved theoretical predictions for the specific 
difference between the ground-state hyperfine splitting values 
of H- and Li-like Bi ions and for the g factor of Li-like Pb ion. 



The screened radiative coiTection in the presence of an ex- 
ternal potential coiTesponds to the third-order perturbation 
theory terms. Nowadays, several approaches are used for 
derivation of the formal expressions from the first princi- 
ples of QED: the two-time Green-function method fllTT . the 
covariant-evolution-operator method fis"], and the line profile 
approach 1 19]. Here, we employ the two-time Green-function 
method. To simplify the derivation of formal expressions, we 
specify the formalism regarding the closed shell electrons as 
belonging to a redefined vacuum. It implies a modification 
of the iO-prescription in the electron propagator incorporat- 
ing the closed-shell electrons. The coiTesponding shift of the 
Fermi-level does not affect the hyperfine structure and the g 
factor. In this way we have to consider all two-loop diagrams 
for the valence electron in the presence of magnetic pertur- 
bation, this is 27 nonequivalent diagrams. These diagrams 
merge the second-order interelectronic-interaction correction, 
the two-loop, and the screened one-loop radiative corrections. 
The generic types of the resulting screened self-energy dia- 
grams are depicted in Fig. [T] The radiative correction to the 
hyperfine splitting Xiad can be written as the sum Xiad = 
2;qed + a^SQED, where xqed corresponds to the one-electron 
QED correction, and xsqed stands for the screened radiative 
coiTection. The latter can be divided into self-energy (SE) and 



vacuum-polarization (VP) parts, xsqed = ^^sqed + ^sqed- 
The screened SE coiTection can be distinguished according to 
so-called irreducible (irr) and reducible (red) parts. It appears 
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FIG. 1: Feynman diagrams representing ttie screened self-energy 
correction in the presence of an external potential. The wavy line in- 
dicates the photon propagator and the double line indicates the elec- 
tron propagators in the Coulomb field. The dashed line terminated 
with the triangle denotes the interaction with magnetic field. 



as the sum of the following terms: 
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Here, a and b denotes the valence and core electron states, 
respectively, the sum over b runs over all closed-shell states, 
P is the permutation operator, giving rise to the sign (—1)^ 
of the permutation, and the notation (a b) stands for the 
contribution with interchanged labels a and b. The SE oper- 
ator E(£), the interelectronic-interaction operator /(w), and 
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TABLE I: Individual contributions to tiie screened SE correc- 
tion a^sQED for the ground-state hyperfine structure of the Li-like 

209gj80+ 



Contr. 


Feynman 


Coulomb 


Contr. 


Feynman 


Coulomb 


A, irr 


0.001544 


0.001555 


F, irr 


-0.000174 


-0.000172 


B, irr 


-0.000380 


-0.000398 


G,red 


-0.001298 


-0.001307 


C, irr 


0.001928 


0.001952 


H,red 


0.000331 


0.000331 


D, irr 


-0.000936 


-0.000945 


I, red 


0.000066 


0.000066 


E, irr 


0.000028 


0.000028 









Total 0.001109 0.001109 



TABLE IL Individual contributions to the specific difference A'i? 
for ^°^Bi in meV. 











A'E' 


Dirac value 




844.829 


876.638 


-31.809 


Interel. inter., ~ 


\IZ 


-29.995 




-29.995 


Interel. inter., ~ 


andh.o. 


0.25(4) 




0.25(4) 


QED 




-5.052 


-5.088 


0.036 


Screened QED 




0.194(6) 




0.194(6) 


Total 








-61.32(4) 



their derivatives (indicated by primes) are defined similar as 
in Ref. lll7ll . u = 1 — zO preserves the proper treatment of 
poles of the electron propagators. The energy difference A is 
defined as A = — £pa and, accordingly, A = e;, — epb in 
terms (a h). Tq is the electronic part of the hyperfine- 
interaction operator and Ga is the multiplicative factor de- 
pending on the quantum numbers of the valence electron (see 
for details Ref. (l5W ). The wavefunctions |^), \(), and |^'), 
IC) are defined as follows 



\n){n\To\a) 



b\PaPb 



\n){nb\I{A)\PaPb) 



(10) 



(11) 



and 1^^) = gf^ICa), KlPaPb) = di:\^b\PaPb/ 

The expressions compiled in Eqs. ([T]i-([3]l, and (|5j-(l8]l con- 
tain ultraviolet divergences. We separate out the divergent 
zero- and one-potential terms in Eqs. ([T]i, (|2]i, (|5]l, (|7]) and 
zero-potential terms in Eqs. (O, (|6]l, (O and evaluate these 
terms in the momentum space, where the divergences can be 
removed analytically (see, e.g., Ref. |(20l]). The remaining 
many-potential terms are ultraviolet finite and calculated in 
coordinate space. The infrared divergences which occur in 
the terms of the Eqs. (|3]l, (|6]l, (O, (|9]l are regularized by 
introducing a nonzero photon mass and canceled analytically. 

The numerical evaluation is based on employing the dual- 
kinetic -balance finite basis set method 12111 with basis func- 
tions constructed from B-splines ll22ll . The Fermi model for 
the nuclear charge density and the sphere model for the mag- 
netic moment distribution have been employed. In what fol- 
lows we present our result for the case of Li-like Bi utilizing 
the corresponding values for the nuclear properties: (r^^/-^ — 
5.5211 fm [23], I"" = 9/2-, and n = 4.1106(2)/i7v H. The 
calculations have been performed in Feynman and Coulomb 
gauges for the photon propagator describing the electron- 
electron interaction, thus providing an accurate check of the 
numerical procedure. The obtained results for the screened SE 
correction for the hyperfine splitting of the Li-like Bi are pre- 
sented in TableUin both gauges, respectively. Finally, we have 
calculated the screened SE correction within local screening 
potentials: Kohn-Sham 0.0012 and core-Hartree 0.0013. The 



results are in reasonable agreement with the rigorous evalua- 
tion ^ 0.00111. 

We have also calculated the screened VP correction in the 
presence of a magnetic field employing the Uehling approx- 
imation for the VP loop. The results have been checked uti- 
lizing the Feynman and Coulomb gauges for the photon prop- 
agator mediating the interelectronic interaction. The electric- 
loop part of the screened Wichmann-Kroll (WK) contribution 
has been calculated by means of the approximate formulas for 
the WK potential from Ref. 12511 . As concerns the screened 
WK magnetic-loop part we have employed the hydrogenic 2s 
value from Ref. [26.1 . assuming that it enters with the same 
screening ratio as the Uehling terms. Accordingly, our value 
for the screened VP correction is a;gQED ~ 

-0.00054(2). 

Finally, the total value for the screened QED correction to 
the ground-state hyperfine structure in Li-like Bi results as 
a^SQED = 0.00057(2). 

Probing the influence of QED effects on the hyperfine split- 
ting of highly charged ions is impeded by the uncertainty of 
the nuclear magnetization distribution correction [the Bohr- 
Weisskopf (BW) effect]. In this context, it was proposed to 
consider a specific difference of the ground state hyperfine 
splitting in H-like and Li-like ions [10]: A'E = AE^'^"') - 
CA£;(i"), where A£;(i") and AE^^^^) are the hyperfine split- 
tings of H- and Li-like ions, respectively, the parameter ^ is 
chosen to cancel the BW correction. In this specific differ- 
ence the nuclear corrections almost vanish completely. We 
have recalculated the energy shift A'E employing the most 
accurate result obtained for the screened QED correction. The 
interelectronic-interaction corrections have been evaluated to 
first-order in 1/Z within the QED perturbation theory and to 
higher-orders within the large-scale configuration-interaction 
Dirac-Fock-Sturm method. Extracting numerically the con- 
tribution of the BW effect in different terms we found that 
the cancellation appears with ^ chosen to be ^ = 0.16886 
for the case of Bi. In Table we present obtained result 
for the specific difference between the hyperfine structure 
values of H- and Li-like Bi ions. The Dirac value incor- 
porates also the finite nuclear size correction. The nuclear- 
polarization correction to the Is hyperfine splitting calculated 
in Ref. (zf] yields S.AE^^p'' = 0.009 meV However, since all 
nuclear corrections have the similar scaling dependence upon 
the principal quantum numbers, we expect the strong cancel- 
lation between Is and 2s nuclear-polarization corrections in 
the specific difference. The same is valid for the second- 
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TABLE III: Individual contributions to tiie ground-state g factor of 
Li-like 2°8pb^9+. 



Dirac value (point nucleus) 


1.932 002 904 


Finite nuclear size 


r\ r\r\r\ mo c o / i o \ 

0.00007s 58(13) 


Interel. inter. 


0.002 140 7(27) 


QED, ~ Q 


0.0024117(1) 


QED, ~ 


- 0.000003 6(5) 


Screened QED 


- 0.0000016(1) 


Nuclear recoil 


0.000000 25(35) 


Nuclear polarization 


- 0.000000 04(2) 


Total 


1.936 628 9(28) 



interelectronic-interaction correction will be the next step to- 
wards the unprecedented accuracy for the stringent test of the 
bound-state QED in the presence of magnetic fields. 
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order one-electron QED contributions. Comparing with the 
results for the specific difference IS!E presented in Ref. [il 
we have increased the accuracy for the screened QED part and 
performed more elaborate calculations for the higher-order 
interelectronic-interaction correction. Further rigorous eval- 
uation of the higher-order electron-electron interaction cor- 
rections will provide a test of bound-state QED at strongest 
electric and magnetic fields. 

Similar calculations have been performed for the g factor of 
Li-like heavy ions. Here, we present our results for the case 
of 208p|^79+ .^jfjj fjjg following value for the nuclear charge 
radius (r2}i/2 ^ 5.501O fm |23]. The rigorous evaluation of 
the screened SE correction gives AggqgQ = —3. 

1(1) X 10-6. 

The previous value obtained with local screening potentials 
was -3.5(1.2) X lO"*" iH- Thus, the uncertainty of the 
screened SE correction has been reduced by an order of mag- 
nitude. The screened VP contribution has been calculated 
within the Uehling approximation. As to the WK part, we 
have employed the approximate formulas for the electric-loop 
potential [25], while the magnetic-loop value has been taken 
from Ref. |28], assuming the same screening ratio as for the 
Uehling term. Accordingly, we have obtained Agggg^ = 
1.5 X 10^6. In Table Hn] we have updated the value for the 
g factor of Li-like 208p|^79+ previously reported in Ref. \\% 
employing the result obtained for the screened QED correc- 
tion. Further extensions of these calculations to the g factor of 
B-like heavy ions may serve for an independent determination 
of the fine structure constant from QED at strong fields [12]. 

In summary, we have rigorously calculated the screened 
QED correction to the hyperfine splitting and g factor of 
heavy Li-like ions. We have increased the theoretical accu- 
racy for the specific difference between the hyperfine splitting 
values of H- and Li-like bismuth as well as for the g factor 
of Li-like lead. The rigorous calculation of the higher-order 



[1] I. Klaft et at., Phys. Rev Lett. 73, 2425 (1994). 

[2] J. R. Crespo Lopez-Urrutia et al, Phys. Rev. Lett. 77, 826 

(1996) . 

[3] J. R. Crespo Lopez-Urrutia et al, Phys. Rev A 57, 879 (1998). 

[4] P Seelig etal, Phys. Rev. Lett. 81, 4824 (1998). 

[5] P Beiersdorfer et al, Phys. Rev. A 64, 032506 (2001). 

[6] H. Hiiffner et al, Phys. Rev. Lett. 85, 5308 (2000). 

[7] J. Verdu et al, Phys. Rev. Lett. 92, 093002 (2004). 

[8] D. R A. Winters et al. Can. J. Phys. 85, 403 (2007). 

[9] M. Vogel et al, Eur. Phys. J. Special Topics 163, 113 (2008). 
[10] V. M. Shabaev et al, Phys. Rev Lett. 86, 3959 (2001). 
[11] D. A. Glazov et al, Phys. Rev. A 70, 062104 (2004). 
[12] V. M. Shabaev et al, Phys. Rev Lett. 96, 253002 (2006). 
[13] J. Sapirstein and K. T. Cheng, Phys. Rev. A 63, 032506 (2001). 
[14] D. A. Glazov et al, Phys. Lett. A 357, 330 (2006). 
[15] A. V. Volotka et al, Phys. Rev. A 78, 062507 (2008). 
[16] H. Persson et al, Phys. Rev. Lett. 76, 1433 (1996); S. A. Blun- 
dell, K. T. Cheng, and J. Sapirstein, Phys. Rev. Lett. 78, 4914 

(1997) ; V. M. Shabaev et al, Phys. Rev. A 56, 252 (1997); 
V. A. Yerokhin and V. M. Shabaev, Phys. Rev. A 64, 012506 
(2001). 

[17] V. M. Shabaev, Phys. Rep. 356, 119 (2002). 
[18] I. Lindgren, S. Salomonson, and B. Asen, Phys. Rep. 389, 161 
(2004). 

[19] O. Yu. Andreev et al, Phys. Rep. 455, 135 (2008). 

[20] V. A. Yerokhin and V. M. Shabaev, Phys. Rev. A 60, 800 (1999). 

[21] V. M. Shabaev et al, Phys. Rev. Lett. 93, 130405 (2004). 

[22] I. Sapirstein and W. R. lohnson, J. Phys. B 29, 5213 (1996). 

[23] I. Angeli, At. Data Nucl. Data Tables 87, 185 (2004). 

[24] N. J. Stone, At. Data Nucl. Data Tables 90, 75 (2005). 

[25] A. G. Fainshtein, N. L. Manakov, and A. A. Nekipelov, J. Phys. 

B 23, 559 (1990). 
[26] A. N. Artemyev et al, Phys. Rev A 63, 062504 (2001). 
[27] A. V. Nefiodov, G. Plunien, and G. Soff, Phys. Lett. B 552, 35 

(2003). 

[28] R. N. Lee et al, Phys. Rev. A 71, 052501 (2005). 



